Time dependence of transmission in semiconductor superlattices 
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Time delay in electron propagation through a finite periodic system such as a semiconductor 
superlattice is studied by direct numerical solution of the time-dependent Schrodinger equation. We 
compare systems with and without addition of an anti-reflection coating (ARC). With an ARC, 
the time delay is consistent with propagation at the Bloch velocity of the periodic system, which 
significantly reduces the time delay, in addition to increasing the transmissivity. 
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I. INTRODUCTION 

Electron transport in a layered semiconductor super- 
lattice (SL) can in many cases be treated as a one- 
dimensional finite periodic system. With as few as TV = 
five cells, a well-developed band structure ensues 
Within each allowed band, where the Bloch phase (j) in- 
creases by TT, the transmission shows narrow peaks de- 
termined by N(j){E) = mTT, m = 1, 2, • ■ • iV — 1. Between 
these peaks the transmission touches an envelope of the 
minima which also has a simple description. In the for- 
bidden bands, the Bloch phase (j) acquires an imaginary 
part; as a result the transmission goes rapidly to zero. 
Pacher et al. ^ have used this property to design an 
electron band-pass filter. Further, by adding a quarter- 
wave cell at each end of the periodic array they 
were able to increase the average transmission within the 
band from about 25% to about 75%. 

In this work we will study the time dependence of an 
electron wavepacket passing through a superlattice, com- 
paring the situation with and without an anti-reflection 
coating (ARC) . In a series of papers [|[0J|,^ , some of the 
present authors have shown how to design an ARC which 
gives optimal transmission within a given miniband, by 
adding suitably configured potential cells on each end of 
a periodic array. The design depends on the shape of the 
potential cells making up the periodic array, but not their 
number. An r-cell ARC consists of r distinct potential 
cells on each end ||]. In the simplest case the periodic 
array consists of refiection symmetric cells; the ARC can 
also be made reflection symmetric by using symmetric 
cells and reversing their order at the opposite end. 

In the transfer matrix formalism, the electron wave 
function at fixed position x is represented by two am- 
plitudes, which can be treated as a spinor. The upper 
component corresponds to right-moving and the lower to 
left-moving waves. Passing though an arbitrary potential 
cell is described by the action of a transfer matrix on this 
spinor. For reflection symmetric cells, the transfer ma- 
trix can be represented in the Kard form |§J9,nG|, which 



involves just two real parameters at given energy. One of 
these is the Bloch phase, and the second, the impedance 
parameter fj,, is approximately constant over the mid- 
dle of an allowed band, but diverges at the band edges. 
Adopting the Kard parameterization allows one the ben- 
efit of a powerful analogy to the precession of an electron 
spin in a magnetic field. The magnetic field direction 
has polar angle 9, where tan6'/2 = tanh/^/2. The angle 
of precession is twice the Bloch phase: 20. An electron 
moving to the right corresponds to spin-up along OZ, 
and left-moving waves to spin-down. Passing through N 
identical cells therefore amounts to precession by angle 
2iV(/); when this is an integer m multiple of 27r, the final 
electron state will be the same as the initial plane wave, 
(except for a phase,) which gives perfect transmission. 

A Bloch state is one where the wave function on ei- 
ther side of the potential cell differs only by the Bloch 
phase (j). This is analogous to an electron whose spin is 
aligned along the direction of the magnetic field. The 
action of an ARC can therefore be understood as taking 
an initial spin-up state and rotating it to lie along the 
field direction. That is, it converts a right-moving wave 
into a Bloch state of the periodic array. Passing through 
any number of cells only adds an overall phase N(j), and 
the downstream ARC reverses the alignment back to the 
spin-up state when the resonant condition holds. 

Without an ARC, electrons which are transmitted do 
so via narrow resonances as mentioned above, so one ex- 
pects a significant time delay. Resonant states and their 
characteristic time dependence [ pdp^ are key parts of 
nuclear physics. Perhaps the definitive discussion of the 
time evolution of wave packets in the vicinity of resonant 
states is that of Rosenfeld When the ratio T/Er of 
resonance width to energy is small, and the wave packet 
is narrow in energy compared to F, there is a broad range 
of intermediate times over which the system exhibits ex- 
ponential time dependence, and therefore significant time 
delay compared to free propagation. The work of More 
and Gerjuoy should also be mentioned. While the 

3D system differs in important ways from ID, the same 
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conclusion holds. The SL is an interesting system to 
study precisely because placing A'^ — 1 resonances within 
an allowed band produces resonances with widths the or- 
der of a few meV, satisfying one of Rosenfeld's criteria. 

With an ARC, propagation via a Bloch state should be 
much faster. To quantify this effect we have performed 
numerical solutions of the time-dependent Schrodinger 
equation (TDSE). The subject of time delay in passing 
though a barrier is a controversial one (see the reviews 
by Hauge and St0vneng Leavens and Aers |T7|,and 
de Carvalho and Nussenzveig ||l8| for example) but it 
is not our purpose to debate the merits or demerits of 
the various definitions that have been used (phase time, 
dwell time, Larmor time, etc.) This controversy con- 
tinues in the context of wanting to define a single time 
to characterize the process. Our approach is to solve 
the time dependent Schrodinger equation directly. This 
provides a direct means of comparing the two situations 
(with/without an ARC), while providing much more in- 
formation than a single number. 
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FIG. 1. Five cell array with a single-cell ARC, comparing 
the original square barrier cells to gaussian barrier cells. 

There have been several papers which discuss prop- 
agation in a superlattice without an ARC. Stamp and 
Mcintosh made a careful study of the two-barrier 
case, which exhibits narrow quasi-bound state (QB) res- 
onances. Their emphasis was on the role of the width (in 
energy) of the incident wave packet on exciting the quasi- 
bound states. They expanded the initial wave packet in 
stationary states of the scattering problem and propa- 
gated the solution in time by quadratures. This method 
is feasible for square barrier/well type potentials where 
the solutions are analytic. Pereyra used a similar 
approach for the superlattice problem. Bouchard and 
Luban considered electrons trapped in an infinite SL, 
in presence of an external electric field. Their empha- 
sis was on the Wannier-Stark ladder of states localized 
by the applied field, and on finding Bloch oscillations. 
Their numerical method is similar to ours in us- 

ing the implicit method to propagate the state forward 
in time. However, they confined the system in a box and 
had to limit the total time so as to avoid refiection from 



the walls. In our work, transparent boundary conditions 
remove such restrictions. The work closest to ours is by 
Pacher and Gornik |4|J^, who considered the effect of an 
ARC on transmission through a superlattice, following 
Pereyra for the time evolution. 

The numerical method we use was pioneered by Gold- 
berg et al. ||2^, and greatly improved by Moyer |25| who 
implemented the Numerov algorithm for the spatial vari- 
ation and added transparent boundary conditions. The 
time-dependence is handled by the lowest-order Crank- 
Nicolson (implicit) method. The wave equation was fol- 
lowed for typically 15 ps, on a region (the "system") 
XL < X < xr of width 3 to 8/im. Transparent boundary 
conditions were especially important in obtaining our re- 
sults; they were applied at both ends of the system, so 
that as the wave function reaches those limits, it will pro- 
ceed outwards without refiection. By integrating the out- 
going fiux at those boundaries, the reflection and trans- 
mission probabilities are accumulated. 

The initial state is a gaussian wave packet sufficiently 
wide in real space to correspond to a small uncertainty 
ctb in energy. For oe/E ~ 3% the root mean square 
(rms) width ax is 140 nm; for <jeIE ~ 0.4% it is 1/zm. 
The main requirement on the system width xr — xl is to 
accommodate the initial wave packet, without overlap- 
ping the potential array. Narrower (in energy) packets 
would require a wider system and more computer time 
for the simulations. 



II. THE CALCULATIONS 

A. Wave packet transmission 

The potential array is located on a < x < b near the 
origin, very close to xji. For this work we adapted a case 
previously studied. It corresponds to a GaAs/AlGaAs 
superlattice of Pacher et al. , but with five barriers rather 
than the six of their device. Their modelling gave the 
barriers a height of approximately 290 meV, and width 
2.54 nm, separated by wells of width 6.50 nm. For con- 
venience we used a constant effective mass m* = 0.071 
in all layers. Further, to simplify the numerical work, we 
replaced Pacher 's square barriers by equivalent gaussian 
shaped barriers. 



d/2 <x < d/2. 



(1) 



The gaussian barrier has height Vq = 270.084 meV and 
a width parameter w = 1.02 nm. The full cell width 
was set at d = 9.4 nm, only 4% wider than the square 
barrier cell. These parameters were fitted to make the 
single-cell transfer matrices equivalent within the lowest 
allowed band, which runs from 50 < < 74 meV. The 
original and the gaussian potentials are shown in Fig. |^. 

By equivalent, we mean that their scattering proper- 
ties are accurately the same, across the allowed band of 
interest. Since the transfer matrix for a symmetric cell 
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depends on only two Kard parameters, we need only fit 
these two, and the results are shown in Fig. ||. The 
cos (j){E) are virtually identical for the original and the 
gaussian potential cell, as seen in Fig. ||(a). The third 
(dotted) line is the cos <j)A of the optimal single-cell ARC. 
That gaussian has height Va — 135.64 meV and width 
parameter wa = 0.98 nm, with a total ARC cell width 
of 7.62 nm. Also was very close for both cells: see 

Fig. ^(b). (The lower line is i^a for the ARC layer. For 
a single-cell ARC the prescription is fiA = 0.5/x at the 
centre of the allowed band.) Having gaussian shape cells 
allowed us to use the Numerov method without having 
to worry about points of discontinuity of the potential. 
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tion (TDSE) we chose gaussian wave packets. The initial 
state 



V'(x,t 0) 
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(2^ai)i/4 
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is normalized to one particle. For a potential array sit- 
ting near the origin, the initial position xq has to be 
taken sufficiently negative so that the wave packet is well 
away from the periodic potential at time zero. For the 
more time consuming runs that followed the transmis- 
sion across the entire band, a width ax = 1000 nm was 
used. When producing videos of the scattering process, 
a less demanding ax ~ 140 nm was chosen. Most of the 
work was done with xl ~ —7500 nm to the left of the 
potential array, and a gaussian wave packet centered at 
the mid-point of that range. For the widest wave packets 
employed, the amplitude of the wave at a;^ and was 
therefore 0.0297 for ax = 1/^m, and was truncated to 
zero. The truncation introduces some high momentum 
components which have a small effect on time develop- 
ment, which is not visible in any of the graphs of this 
paper. Had we doubled the distances to 15 /zm, the trun- 
cation would have been at 8 -10^^ in amplitude and not 
at all discernible. The fractional standard deviation in 
energy can be written as 



aE/E 



4.24nm /60 meV 



E 



(3) 
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b) 



FIG. 2. Comparison of (a) Blocli pliases of a square barrier 
cell (dash line) and our gaussian cell (solid line); also shown 
is cos <^A of a single-cell ARC (dotted line) . (b) same for 
impedance parameters fj,, /ia- 

In Fig. H we show the transmission profiles of the 5-cell 
SL with no ARC (dotted line), a single-cell ARC (dashed 
line) and a two-cell ARC (solid line), for the gaussian 
cells array (computed using the usual time-independent 
methods). With the two-cell ARC, the high transmission 
region runs from 53 to 68 meV, and the satellite peaks 
are pushed closer to the band edges. This agrees well 
with the square barrier calculations in 

For solution of the time-dependent Schrodinger equa- 
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FIG. 3. Transmission for 5-cell array, bare and with sin- 
gle-cell and double-cell ARC (time-independent calculation). 

In Fig. ^ we show transmission for the gaussian ar- 
ray computed using the TDSE. The right-moving flux at 
xr = X is 



2ime'm* ^ ' ' ^j^, 
and the integrated transmission is 



F (2^'*) — 7Z i^ix^tJl (4) 



T{Ep,t)= I j{x,t)dt with 
Jo 

Tas{Ep) = lim T{Ep,t) . 

t — '■OO 



(5) 
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The label Ep is mean energy of the wave packet. 
Corresponding expressions for the reflection probabil- 
ity R{Ep,t) hold with the left-moving flux monitored 
at XL = X. For packets narrow in energy T^s 
should approach the transmission probability of the time- 
independent solutions, and this is seen to occur. The 
main difference with respect to Fig. || is that the peaks 
are smeared by the finite width of the wave packet (ap- 
proximately 0.5 meV). Even so, they are so narrow that 
the tops are ragged due to the finite steps in Ep of the 
wave packets used to generate the figure. (With finer 
steps and narrower (in energy) wave packets, convergence 
has been verified.) 




Also in Fig. lines below the peaks show their build- 
up over time T{E,t)] this proceeds uniformly, with some 
bias to faster development at higher energies. This bias 
can be understood as a simple velocity effect: the higher 
energy wave packets travel faster and reach the right wall 
xr sooner than the slower moving wave packets near the 
lower band edge. We note an apparent discrepancy with 
Fig. 4 of Romo |2^: at 0.6 ps he shows transmission in 
the troughs exceeding the asymptotic {t > 60 ps) trans- 
mission. However, as pointed out by an astuute reader, 
Romo's calculation is based on a very different initial 
state than ours |^,^. At time t = {) his initial state 
is a uniform standing wave confined by a mirror to the 
left of the array at a; < a. The mirror is removed at 
t = 0, and after an infinite time a new equilibrium state 
is achieved which is the stationary scattering state for a 
wave incident from the left. As the wave front advances 
to the right, the leading edge becomes smeared out, and 
oscillations develop behind, a process which Moshinsky 
prf called diffraction in time. Ultimately the limit as 
t — > oo of |7/;(6, t)p should approach the time- indepen- 
dent transmission probability T{E). The curves at inter- 
mediate times, at a fixed position, show the passing of 
the wave front and the characteristic oscillations behind 
it. They do not represent the accumulation of T(Ep,t) 
as in our calculation. Indeed, we can calculate using a 



long square wave packet normalized to one particle, and 
we find similar oscillatory behaviour in the build-up of 
both \tp{x,t)\'^ and the flux, which is largely absent in 
the time-integrated flux. 

Fig. 1^ shows similar results for the same array plus a 
single-cell ARC. Again the asymptotic transmission is in 
good agreement with a time- independent calculation, but 
smeared by the finite width of the wave packet. Also, the 
build-up of the transmission profile proceeds smoothly 
with the above-noted velocity bias. 
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FIG. 5. Transmission for gaussian array plus single-cell 
ARC, showing development over the same time intervals. 
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FIG. 6. Energy content of the transmitted, reflected, initial 
and final wave packets near the 51 meV resonance. 

Near each band edge there is a satellite peak which 
resembles the narrow resonances of the periodic array. 
To study the character of these lines we scanned over 
them using spatially narrow wave packets. We enlarged 
the system to the right, so that there was room for both 
the reflected and transmitted wave packets to become 
well separated from the array. We then Fourier analysed 
the reflected and transmitted wave packets individually. 
The Fourier transforms (moduli squared) R{E) and T{E) 
are plotted against energy in Fig. || and the results are 
quite revealing. The T packet is very narrow compared 
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to the incident wave, and centered just below 51 meV. 
The R packet in contrast has a node at this energy, and 
corresponds to the energies on either side. This shows 
that the transmission does proceed primarily through a 
narrow resonance. The transform (mod squared) of the 
complete packet after scattering, denoted "sum R+T" 
sits right on top of the curve for the initial wave packet, 
which is a testament to the unitarity of the numerical 
work. 

The time-development of the reflection and transmis- 
sion probabilities is shown in Fig. ^ for a narrow (in 
energy) wave packet centered at i?p = 51 ± 0.47 meV. 
T{t) and R{t) (see eq. |) are accumulated as the waves 
exit from the right/left boundaries xr,xl- T{t) begins 
to rise at 2 ps; this corresponds to the time taken for part 
of the initial packet to be transmitted and reach x^. The 
longer time scale for R{t) is due to the wider space on the 
left which accommodated the initial wave packet. One 
sees two steps in the R{t) curve, one corresponding to 
direct reflection from the leading edge of the array, and 
the second one from entrapment before reflection. The 
sum R{t) + T{t) ~ 1 — P{t) is the complement of the 
probability P{t) remaining within the system (x^jX^): 



P{t) = / \i^{x,t)\^dx . 

J XL 



(6) 



(A 60 meV electron with effective mass 0.071 travels at 
545 nm/ps; the time offsets are consistent with this.) 
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FIG. 7. Accumulated transmission (T), reflection (R) and 
probability of remaining in the system (P) versus time, for a 
narrow wave packet at iJp = 51 ± 0.47 meV. 



B. Quasi-Bound states 

At times near 6 ps T{t) is approaching its asymptotic 
value Tag. Looking at the wave which is still trapped 
inside the array, one sees that it has a particular shape 
(shown in Fig. H) and is decaying with mean life Tn = 
0.838(2) ps. This corresponds to a half-width F„ = 0.786 
mcV. It implies a quasi bound state located at zF„/2 
in the lower half complex plane. Without the ARC, the 



position would be slightly displaced, corresponding to the 
shift of the outer peaks in Fig. ||. 

An electron inside the array sees four potential wells 
between the five stronger barriers. The QB states are 
built on the ground state in each well. Through coupling 
across the barriers four QB states are created, associated 
with the four transmission peaks seen in Fig. |^ for ex- 
ample. Fig. H shows that even with the ARC, the first 
of these states persists (as does the last). Rosenfeld |l^] 
discussed in detail the conditions under which the expo- 
nential decay of such states would be observed; the most 
important is that the ratio r„/i?„ << 1, which is well 
satisfied for our superlattice. Stamp and Mcintosh |l^] 
emphasized another criterion, the "interaction time" of 
the wave with the potential. This is basically the dwell 
time inside the potential. If it is too long compared to 
the lifetime of the QB state, the state will be decaying 
continuously as it is being fed, and the probability will 
not be built up to a significant extent. For the state of 
Fig. |8[ Tn is of the same order as the time delay shown 
in Fig. Il^, below. 
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FIG. 8. Decay of the QB state at 50.8 meV. The potential 
array with single-cell ARC (solid line) is shown for orientation. 

The fourth QB state also survives the addition of the 
ARC, and implies a pole at 72 — i0.767 meV. The draw- 
ing would look very similar to Fig, except that the 
wave function has nodes in each barrier. That is because 
the lowest state is nodeless, while the fourth state has 
alternating signs in the four wells. 



III. TIME DELAY 

A. Results from time-dependent calculations 

The primary motivation for this work was to compare 
the time taken to traverse the array plus ARC, as com- 
pared to the simple array. However, in the case with 
ARC, the two satellite peaks appear to be narrow res- 
onances of the same type as for the simple periodic ar- 
ray. Hence it is also of interest to compare transmis- 
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sion through the satelhte peaks with transmission via 
the Bloch states of the broad central maximum. 

Fig. shows that the build-up of the transmission 
peak at 51 meV proceeds linearly in time, over the range 
0.2 < T{t) < 0.6. Similar linear behaviour applies 
throughout the allowed band. In the absence of a po- 
tential a similar curve is found, with some time displace- 
ment. By comparing the two, a reliable time delay (or 
advance) can be deduced which depends little on just 
which point is selected. For larger and smaller times the 
build-up is non-linear, and the converse is true. 

In Fig. ^ (a) we have plotted the ratio of T[t)/Tas 
(solid line) as a function of energy across the whole band. 
The several lines correspond to different times. The in- 
clined straight dotted lines are the same thing in the 
absence of any barriers (free propagation). In the back- 
ground as a chain line, the transmission probability T{E) 
is plotted for reference. Except at the four transmission 
peaks, the solid and dotted lines agree well. Also, the 
rate of build-up is the same as in the absence of the bar- 
riers. That says that in the transmission troughs, almost 
nothing is being transmitted but it does so with no time 
delay. Under the four peaks, the build-up though still 
linear in time, is significantly retarded compared to free 
propagation. The inclination of the dotted lines we in- 
terpret as a velocity effect: the higher energy electrons 
travel faster and arrive at xr sooner. 

To deduce a time delay, we take the difference be- 
tween the solid and dotted lines, averaged over the range 
0.3 < T{t)/Tas < 0.7. This avoids using the non- linear 
portion of the T{t) curve, as already discussed in con- 
nection with Fig. 1^. It would make little difference had 
we simply taken the time delay at ratio 0.5 . The result- 
ing time delay is plotted in Fig. |l^ (a). This procedure 
is consistent with the work of Dumont and Marchioro, 
who argue that for suitably chosen wave packets and by 
measuring the flux of particles at x^, one can define a 
"tunneling time probability distribution" , which is ba- 
sically the derivative of the integrated transmitted flux 
T{t) plotted in Fig. |^, but normalized to its asymptotic 
value Tas- 

Fig. H (b) shows the same ratio T{t)/Tas after adding 
a single-cell ARC. Except at the two satellite peaks, the 
solid and dotted lines are roughly parallel, but not touch- 
ing, showing that there is a delay for propagation via a 
Bloch state. The delay is significantly greater under the 
satellite peaks. The corresponding time delays are plot- 
ted in Fig. |l^ (b). Without an ARC, they go to zero 
(or even negative) in the transmission valleys, but at the 
peaks they range from 0.75 to 2.6 ps. With a single-layer 
ARC the delay is 0.25 ps over the middle of the band, 
jumping to 1.6 and 0.8 ps under the satellite peaks. So, 
for those electrons that are transmitted, adding the ARC 
cuts the time delay in half, while greatly increasing the 
average transmissivity. [Also shown in Fig. |l^ (b) is a 
simple estimate of the time delay for traversing five cells 
at the Bloch velocity, ignoring any delay in the ARC cells. 
It can be seen that this estimate is the right size, and the 



downward trend from left to right also agrees. Because 
the Bloch velocity vanishes at a band edge, the curve di- 
verges there. Calculations with ten and fifteen cells show 
that the satellite peaks move closer to the divergence at 
the band edge. The delay at 51 meV, 1.6 ps, is about 
double the lifetime of the quasibound state, but at 70 
meV the two are about equal. Similar conclusions hold 
when a two- cell ARC is added, but for brevity we do not 
include them here. 



B. Comparison with time-independent calculations 

Time delay can also be computed in the time- 
independent formalism. From the vast literature on this 
subject, we base our discussion on relatively recent work 
of Nussenzveig [^,|l^ . He argues in favour of the "av- 
erage dwell time" as the most appropriate measure of 
the time taken to pass a barrier, and lists its advantages. 
It is closely related to the "phase time" or transmission 
group delay hdrj/dE |31 1 originally introduced by Wigner 
and Eisenbud. Here fj is the phase of the complex trans- 
mission amplitude, closely related to the transfer matrix 
element by Mn — 1/t, as we now explain. 

In the transfer matrix method, our wave functions 
are defined with a different phase than generally used: 
namely the phase is set to zero on each side {x = a, b) of 
the potential array, rather than at the origin. Adopting 
Nussenzveig's notation for the asymptotic wave function 
to left and right, we compare our ip{x) with the usual 
convention ipix) as follows: 



V'(x) 



ikx 



+ re-*'=(^-2'^) ; te 



ik{x — 'w)-\ 



(7) 



where w = 6 — a is the total width of the potential. It 
follows that the phase rj of our transmission amplitude is 
related to the usual one hy i) = rj — kw. 
Then the phase time delay is 



df] 

Tph — n-r^ 



dE 



drj 
dE 



w 



dk 
dE 



(8) 



where rj is the phase of our transmission amplitude and 
w is the width of the superlattice (not the entire system) 
under consideration. Since dE/dk is the velocity Vf of a 
free particle, the contribution w/vf represents the time 
taken to cross the superlattice assuming zero potential. 
In view of our earlier estimate for the free velocity, it is of 
the order of 0.1 ps. Thus, dfj/dE gives phase time delay, 
while our drj/dE gives phase time. 

For the case of unilateral incidence on a symmetric 
potential, which we use, Nussenzweig arrives at his eq. 
(20) for the mean dwell time delay as a spectral average 
over the phase time, plus an oscillatory term. In our 
notation, his result can be written as in eq. 
appendix. 
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In Fig. we show the time delay extracted from our 
time-dependent calculations, both before (a) and after 
(b) a single cell ARC has been added. In both panels, 
the lower dashed line would be the phase time delay, 
assuming that in the periodic array the wave propagates 
at the Bloch velocity of the infinite periodic system: 



VBl 



— o?sin0 dE 
h d cos ( 



d/TBi , 



(9) 



where d is the cell size. Within an allowed band, cos cj) is 
generally quite linear in E. 

In Fig. [lO|(b) the dash-dot line includes an estimate of 
the time delay for passing through the ARC layers as well 
as the five central cells (it is a 10% effect). This Bloch 
time-delay agrees quite well with the result of our time- 
dependent calculations, especially in its general trend. 
This confirms our understanding of the action of the ARC 
layers: they convert the incident plane wave state into a 
Bloch state of the periodic system. The divergence at 
the band edge arises because the Bloch velocity vanishes 
there. 
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FIG. 9. Ratio T/Tas{E) with (solid line) and without po- 
tential (dotted line) for (a) five barrier periodic system and 
(b) with single-layer ARC added. The contours are plotted at 
intervals t = 2 (1.5) 14 ps. Also shown (dotted) is the trans- 
mission T{E), for reference. 

In Fig. |l^(a) the dash-dot line is the locus of phase 
time delay at transmission maxima. While maximal time 



delays do not occur exactly at the transmission maxima, 
they do lie close together for this type of potential cell, 
and the difference would only be visible on a magnified 
drawing. The locus passes about 10% above the peaks, 
which we ascribe to two effects: (i) the finite steps in 
mean wave packet energy in our calculations, which may 
miss the top, and (ii) the finite energy width of the wave 
packet which smears out the result based on plane waves, 
as in eq. |l^. This is seen more clearly in Fig. |l^. Inter- 
estingly, the time delay without an ARC oscillates around 
the Bloch time delay in Fig. |l^ (a); this is not partic- 
ularly obvious from the present case of five cells, but if 
one compares systems with ten to fifteen cells it is quite 
striking. 
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FIG. 10. Time delay (solid line) for electrons transmitted 
through (a) the five barrier periodic system and also (b) with 
single-layer ARC added. Also shown is the transmission T{E) 
(dotted line) for reference. 

The locus of time delay at transmission maxima is de- 
rived as follows: For a finite periodic system, we know 
that (I 

^ ^ ■sinA^0-sin(A/'-l)(/)) . (10) 



In 
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We write t = |t|e^'' and tN = |tAr|e^''" for N cells. Thus 
77 jv can be expressed in terms of rj and the Bloch phase 
(f). Using eq. (||) for the time delay we can similarly ex- 
press T in terms of the single cell parameters. To look 
for maxima, we set the derivative with respect to energy 
to be zero, and sin N(j) = 0, because transmission max- 
ima occur when N(j) is an integer multiple of tt (see first 
paragraph of this paper.) The result is 



.dri_ 



'N 



dE 



Nd sin 77 



\t\vBi sini 



= Ntbi cosh ji 



(11) 



where tbi is the time taken to cross one cell at the Bloch 
velocity, and fi is the impedance parameter. The plot 
of cosh /I in Fig. ^ explains the shape and height of the 
curve immediately. 

From eq. |l^ we subtract the last term of eq. (||), 
leading to an expression for the locus of time delay at 
transmission maxima; 



Tloc = Ntbi cosh^ - 



Nd 



Vfr 



(12) 



Results using eq. (^ are in excellent agreement with 
those of the time-dependent calculation, when allowance 
is made for averaging over energy in the neighbourhood of 
the sharp resonances. There the finite width (in energy) 
of our incident wave packet mainly reduces the height of 
the peaks of the time-independent result, by 5 to 10%, 
acording to eq. (p^). An example is shown in Fig. |l^. 
In this drawing, we should have averaged the phase time 
over the spectral content of the initial state. We did 
not, but because our wave packets are narrow in energy, 
convolution only reduces the heights of the narrow peaks, 
which can be estimated from 



9iE) 



f{E) 



l + {2aE/TYf{E) 



(13) 



This assumes that the unsmeared function J{E) is of 
Breit-Wigner form with width F, and is wide compared 
to the width in energy ~ 0.4 meV of the wave packet. 
Eq. (0|) agrees well with the reduction seen in Figs. [lO| 
and nir both for the time delay and transmission peaks. 
Mostly hidden below the solid line in Fig. |ll| is our time- 
delay from Fig.|l^ (a). One can see that except just below 
the peaks, the agreement is excellent. 

Neither have we included the oscillatory term in the 
mean dwell time delay, which arises from interference be- 
tween the incident wave and the reflected wave pl| , p2| . 
In our calculations there are indeed spectacular interfer- 
ence effects within and to the left of the potential array, 
as the reflected wave is generated. But after the reflected 
wave is well separated from the potential and reaches the 
counter at xl, these oscillations have disappeared, and do 
not show up in the integrated flux. Neglect of this term, 
in the present calculation, is justified in the appendix. 

Pacher and Gornik [^|J|] have also computed tunneling 
times using a similar formula following Pereyra pCl] , with 



similar results. They did not perform time-dependent 
calculations to establish the validity of the result. Fur- 
thermore, in our reading, Pereyra's derivation pl| ] as- 
sumes that the reflection amplitude has a fixed modu- 
lus, and only the phase is varying, which is obviously 
questionable in the case of sharp resonances. After this 
work was submitted, Pacher et al. |Q have also discussed 
phase time delay at the transmission maxima, deriving 
the result in eq. 
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and many others. 
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FIG. 11. Time Delay (solid line) computed from phase of 
tisi, and (dashed line) from the time-dependent calculation, 
just visible below the peaks. The dash-dot line is the locus, 
eq. |l2| Also shown for reference is the transmission (dotted 
line) . 



IV. CONCLUSION 

Time dependence of scattering from a finite periodic 
potential array was studied by direct numerical solution 
of the Schrodinger equation. It was verified that the nar- 
row transmission peaks are associated with quasi-bound 
states of the array. These resonances entail time delays of 
order 1 to 2 ps, while in the transmission minima the de- 
lay vanishes. Upon addition of an anti-reflection coating, 
the broad central transmission maximum corresponds to 
transmission via Bloch states, with a time delay of or- 
der 0.2 to 0.3 ps, as seen in Fig. |l^ (b). The satellite 
peaks near the band edge continue to proceed through 
QB states, but even their time delay is cut roughly in 
half, as compared to the bare periodic array. 

Our incident wave packets had widths in energy of or- 
der 0.4 meV; this was feasible due to the application of 
transparent boundary conditions at the edges of the re- 
gion considered. A further improvement in the method 
of solution is possible, by going from first-order to third- 
order Crank-Nicolson integration for the time depen- 
dence p^ , ^ . This would match the truncation error 
of the Numerov method, while greatly speeding up cal- 
culations. Systems under bias of an applied electric field 
pa can also be handled by this method. Extension to 
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two-dimensional systems is also under consideration. 
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APPENDIX A: SPECTRAL AVERAGE OF THE 
DWELL TIME DELAY 

For our wave packet eq. ^, the Fourier transform is 

Ml) = [8^fTi]i/4e-*9=^''e-(«-'^')'"' , (Al) 



normalized according to 



\Mq)\ 



2tt 



(A2) 



The spectral weight function is therefore A(q) = 
■(/'o('z)|^/(27r). Nussenzweig's eq. (20) for the time de- 
lay by a symmetric potential, with unilateral incidence, 
in our notation becomes 



< Atd >f!f""" = 

dfj 



h\riq)\ 



dE 2E{q) 



cos(2(72;i — f]) 



dq 
2^ 



(A3) 



The terms in square brackets are the dwell time delay 
in the mono-energetic case. Razavy ||3^ for example de- 
rived them by following the method of Smith , albeit 
with some typos in his eq. (18.19). (One has to note that 
for a symmetric potential his two phase shifts are related 
by 77 = 5 + ^/2.) 

As stated earlier, we used a wave packet with a width 
ax of order 1000 nm. The width in energy is of order 
0.4 meV, which is small compared both to the width of 
the allowed band, 25 meV, and even of the resonances of 
a finite periodic potential array, which are in the range 
of 2 to 5 meV depending on the position in the band. 
(The narrowest states are those crowded against the band 
edge.) It is reasonable to treat A{q) as narrow compared 
to the width of the peaks in transmission. 

The first term, drj/dE, varies on the scale of the band- 
width divided by TV, the number of cells, about 6 meV 
in our calculation. The slope of 77 is steepest at a reso- 
nance, where 77 — mn, and is minimal at the mid-points 
between resonances. The phase shift 77 varies smoothly 



compared to the wave packet, and that leads to the con- 
clusion stated in eq. y^: the spectral average mainly 
reduces the height of each transmission peak, leaving its 
position and width unchanged. 

The second, oscillatory term, (OT), includes a factor 
cos(2fca::L — vi^))- Since our xl ^ —7500 nm., this is 
a very high frequency oscillation, given that we set the 
boundary xl far to the left of the potential. The re- 
flection amplitude \r{k)\ varies on the same scale as the 
phase shift. In doing the spectral average it is reason- 
able to treat the small prefactor h\r{k)\/ E{k) ~ 10 fs 
as slowly varying in comparison to the cosine. The mean 
value of the OT can be estimated by the method of steep- 
est descents, leading to 

< OT >^ cos(2fca:i - ry(fc)) 

J e^p[{i{q - k){2xL) - (q - kf2al]dq 

h\r{k)\ 



2E{k) 



cos{2kx L - ri{k)) cxp[-(a;L/crx)72] • (A4) 



By taking xl sufficiently far to the left, the spectral av- 
erage can be made as small as we please. Mostly we used 
s^l/cx = 7.5, so the exponential factor is e~^^, multiply- 
ing a term which is already small. 
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